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Abstract 

This paper provides a comprehensive discussion of neutrahno dark matter within 
classes of extended supersymmetric models referred to as the USSM containing 
one additional SM singlet Higgs plus an extra Z', together with their superpart- 
ners the singlino and bino'. These extra states of the USSM can significantly 
modify the nature and properties of neutralino dark matter relative to that of the 
minimal (or even next-to-minimal) supersymmetric standard models. We derive 
the Feynman rules for the USSM and calculate the dark matter relic abundance 
and direct detection rates for elastic scattering in the USSM for interesting regions 
of parameter space where the largest differences are expected. 



1 Introduction 



One of the benefits of weak scale supersymmetry (SUSY) with conserved R parity is 
that the lightest supersymmetric particle (LSP) is absolutely stable, and provides a 
weakly interacting massive particle (WIMP) candidate [1, 2] capable of accounting for 
the observed cold dark matter (CDM) rehc density QcDMh"^ = 0.1099 ± 0.0062 [3]. In 
particular, the lightest neutralino in SUSY models is an excellent candidate, providing 
its mass, composition and interactions are suitably tuned to result in the correct value of 
QcDAdh'^- The minimal supersymmetric standard model (MSSM) has become a widely 
studied paradigm [4]. However the stringent upper bound on the Higgs boson mass 
in the MSSM combined with its experimental lower bound from LEP has led to some 
tension in the electroweak symmetry breaking sector, roughly characterized by a fine- 
tuning of parameters at the percent level [5]. While the experimental elusiveness of the 
Higgs boson may cast some doubt on the MSSM, there are a host of non-minimal SUSY 
models which predict a heavier and/or more weakly coupled Higgs boson [4]. 

A further reason to move beyond the minimal case is the so-called mu problem of 
the MSSM [6]. The MSSM contains a bilinear mass term that couples the two Higgs 
doublets with a dimensionful coupling fi. This term is SUSY preserving, and as such 
only has two natural values, n = and /i = Mpi (unless special forms of the Kahler 
metric are assumed). However experimental data and the stability of the Higgs mass 
requires that // be of the order of the SUSY breaking scale. In non-minimal SUSY models 
the mu problem is solved by setting /x = and including an additional superfield S, a 
singlet under the Standard Model (SM) gauge group, which couples to the Higgs doublet 
superfields Hi, H2 according to XSH1H2, where A is a dimensionless coupling constant. 
We shall refer to such models generically as singlet SUSY models. Such a coupling 
replaces the SUSY Higgs/Higgsino mass term 1JLH1H2 of the MSSM. The singlet vacuum 
expectation value (VeV) {S) then dynamically generates a SUSY Higgs/Higgsino mass 
near the weak scale as required. This results in an increased Higgs boson mass upper 
bound depending on the value of A, and hence a welcome reduction in electroweak fine 
tuning in addition to solving the /x problem of the MSSM [7] . 

However, although an extra singlet superfield S seems like a minor modification to 
the MSSM, which does no harm to either gauge coupling unification or neutralino dark 
matter, there are further costs involved in this scenario since the introduction of the 
singlet superfield S leads to an additional accidental global U{l)x (Peccei-Quinn (PQ) 
[8]) symmetry which will result in a weak scale massless axion when it is spontaneously 
broken by {S) [9] . Since such an axion has not been observed experimentally, it must be 
removed somehow. This can be done in several ways resulting in different non-minimal 
SUSY models, each involving additional fields and/or parameters. For example, the clas- 
sic solution to this problem is to introduce a singlet term 5*^, as in the next-to- minimal 
supersymmetric standard model (NMSSM) [10], which reduces the PQ symmetry to the 
discrete symmetry Z^. The subsequent breaking of a discrete symmetry at the weak 
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scale can lead to cosmological domain walls which would overdose the Universe. This 
can be avoided by breaking the Z3 symmetry explicitly without upsetting the hierar- 
chy problem by non-renormalizable operators that obey a i?-symmetry [11], or by 
removing the 5*'^ term altogether [12]. 

Another solution to the axion problem of singlet models, which we follow, is to 
promote the PQ symmetry to an Abehan V{V)x gauge symmetry [13]. The idea is 
that the extra gauge boson will cat the troublesome axion via the Higgs mechanism 
resulting in a massive Z' at the TeV scale. The essential additional elements of such 
a scenario then consist of two extra superfields relative to those of the MSSM, namely 
the singlet superfield S and the U{\)x gauge superfield B' . The scenario involving 
only the MSSM superfields plus these two additional superfields, may be considered as 
a phcnomcnological model in its own right which has been referred to as the USSM. 
In the USSM, then, the MSSM particle spectrum is extended by a new CP-even Higgs 
boson S*, a gauge boson Z' and two neutral -inos: a singlino S and a bino' B' while other 
sectors are not enlarged. The presence of new singlino and bino' states greatly modifies 
the phenomenology of the neutralino sector both at coUiders and in cosmology-related 
processes. The collider phenomenology and cosmology of the USSM has been studied 
in [14, 15, 16, 17, 18, 19, 20], which we briefly review as follows. 

The collider phenomenology of the USSM has recently been considered in [18]. The 
neutralino production cross sections and their decay branching fractions depend crucially 
on their masses and composition with respect to the MSSM case. If the new -ino states 
are heavy, their influence on the MSSM neutrahnos is small. In contrast, if the singlino 
mass scale is low, the production rates can be quite different and, since there arc more 
ncutralinos, the decay chains of sparticles can be longer. Moreover, if the mass gaps 
between the MSSM and new -inos are very small, the standard decay modes are almost 
shut and radiative transitions between neutralino states with a soft photon may be 
dominant. In such a case the decay chains can be apparently shorter, a feature which 
is of relevance for the LHC experiments. 

The dark matter phenomenology of the USSM was first studied in [15, 16], and more 
recently in [19]. In [15], the analysis was performed in a scenario with a very light Z' and 
considered the case in which the LSP was a very light singlino. This allowed the authors 
to consider the annihilation of dark matter in the early universe to be dominated by 
s-channel Z' processes, allowing an analytic solution to the dark matter relic density to 
be obtained. In the full parameter space of the USSM, this is just one possibility, and 
indeed such a light Z' is heavily disfavored by current data. In [19] the recoil detection 
of the dark matter candidate in the USSM (and other non-minimal SUSY models) was 
considered. 

In this paper, we provide an up to date and comprehensive analysis of neutralino 
dark matter in the USSM^. We provide a complete discussion of the extended gauge, 

^The recent observation of a positron excess by the PAMELA collaboration [21] have caused a flurry 
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neutralino, Higgs squark and slepton sectors in the USSM, and using the LanHEP [24] 
package, derive all the new Feynman rules involving these extended sectors. We first 
provide a complete qualitative discussion of the new annihilation channels relevant for 
the calculation of the cold dark matter relic density for the neutralino LSP in the USSM. 
We also discuss the elastic scattering cross section for the neutralino LSP in the USSM, 
including both spin-independent and spin-dependent parts of the cross sections, relevant 
for the direct dark matter search experiments. We then survey the parameter space of 
the USSM, and discuss quantitatively how the nature and composition of the neutralino 
LSP can be significantly altered compared to that in the MSSM due to the extra singlino 
and bino' states, for different ranges of parameters. The Feynman rules are then imple- 
mented into the micrOMEGAs [25] package in order to calculate the relic density for 
the corresponding regions of parameter space. This provides a full calculation of the 
annihilation channels including co-annihilation and careful treatment of resonances as 
well as accurately calculating the relic density for an arbitrary admixture of states. In 
this way we extend the analysis of USSM neutralino dark matter annihilation beyond 
the specific cases previously studied in the literature. We also perform an equally gen- 
eral calculation of the direct detection cross-sections for USSM dark matter for elastic 
neutralino-nuclci scattering. 

It is worth emphasizing that the USSM is not a complete model, since from its 
definition it does not include the additional superfields at the TeV scale, charged under 
the gauged Abelian symmetry, which are necessarily present in order to cancel the 
fermionic gauge anomalies involving the U{l)x gauge symmetry. For example, a well 
motivated and elegant solution to the problem of anomaly cancelation is to identify the 
Abelian gauge group as a subgroup of Eq and then cancel the anomalies by assuming 
complete 27 dimensional representations of matter down to the TeV scale. With the 
further requirement that the right-handed neutrino carries zero charge under the Abelian 
gauge group (in order to have a high scale see-saw mechanism) this then specifies the 
theory uniquely as the EgSSM [26, 27]. However our working assumption is that the 
additional matter superfields required to cancel anomalies are heavy compared to the 
Z' mass. The USSM considered in this paper may thus be regarded as a low-energy 
truncation of the EgSSM model, with other EgSSM fields assumed heavy, and the charge 
assignments under the extra U{l)x as given in [27] and summarized in Section 2. 

Despite that fact that the USSM must be regarded as a truncation of a complete 
model, it makes sense to study the physics and cosmology in the USSM since it provides 
a simplified setting to learn about crucial features which will be relevant to any complete 
model involving an additional U gauge group and a singlet. For example, as already 
mentioned the neutralino LSP in the USSM may have components of the extra gaugino 

of speculation that the high energy positrons are produced by annihilating dark matter in the galactic 
halo[22]. An alternative explanation is that astrophysical sources coiild account for the positron excess 
- in particular nearby pulsars [23] . It is unclear as yet which of these explanations is correct and as a 
result we do not address the PAMELA results further in this work. 
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B' and singlino in addition to the usual MSSM neutralino states. Naively we might 
expect that the dark matter phenomenology of such regions would be similar to that 
of singlino dark matter in the NMSSM. However this is not the case. The inclusion of 
the bino' state, as well as the lack of a cubic interaction term 5*'^, results in a signif- 
icant change in the phenomenology. Also the neutralino mass spectrum in the USSM 
is very different from that of the NMSSM as the singhno mass is determined indirectly 
by a mini-see-saw mechanism involving the bino' soft mass parameter M[ rather than 
through a diagonal mass term arising from the cubic S^. The lack of a cubic interaction 
term also restricts the annihilation modes of the singlino, making it dominantly reliant 
on annihilations involving non-singlet Higgs bosons and higgsinos. As the USSM has 
a different Higgs spectrum to the NMSSM, notably in the pseudoscalar Higgs sector, 
the Higgs dominated annihilation channels of the USSM singlino are significantly mod- 
ified with respect to the NMSSM singlino. As Higgs exchange diagrams dominate the 
direct detection phenomenology, the difference in the Higgs spectrum and the singlino 
interactions results in significant differences in the direct detection predictions as well. 

The remainder of the paper is organized as follows. In section 2 we shall define the 
Lagrangian of the USSM and discuss the Higgs, Z' , neutrahno and sfermion sectors. 
In sections 3 and 4 we shall give an overview of the important features of the relic 
density calculation and the direct detection calculation, respectively, highlighting the 
main differences to the MSSM. In section 5 we present the results of the full numerical 
calculations for both the relic density and the direct detection cross-section. It will be 
performed in two physically interesting scenarios: (A) with the MSSM higgsino and 
gaugino mass parameters fixed, while the mass of the extra U(l) gaugino taken free (to 
complement the collider phenomenology discussed in Ref. [18]); (B) with GUT-unified 
gaugino masses. Section 6 summarizes and concludes the paper. The mass matrix 
structure of the extended Higgs scalar sector, and a discussion of the Feynman rules 
involving the extended neutralino sector in the USSM are given in a pair of appendices. 

2 The USSM model 

Including the extra U(l)x symmetry, the gauge group of the model is G = SU(3)c x 
SU(2)i X U(l)y X U(l)x with the couplings 513, 512, Qy-i gx, respectively. In addition to the 
MSSM superfields, the model includes a new vector superfield Bx and a new iso-singlet 
Higgs superfield S. The usual MSSM Yukawa terms Wy of the MSSM superpotential 
(i.e. without the term) are augmented by an additional term that couples the iso- 
singlet to the two iso-doublet Higgs fields: 

W^WY + XS{H^Hd). (1) 

The coupling A is dimensionless. Gauge invariance of the superpotential W under U{l)x 
requires the U(l)x charges to satisfy + Q^^ + Qg = and corresponding relations 
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between the U(l)x charges of Higgs and matter fields. In the following we adopt the 
U(l)x charges as in the EgSSM model [27], see Table 1. (For notational convenience 
we will also use Qi = Qh^^ Q2 = Qhu Qs = Qs-) "^^^ effective /j, parameter is 
generated by the vacuum expectation value {S) of the scalar ^-field. 
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Table 1: The U{1)y and U{l)x charges of matter fields in the USSM, where Qf and QJ are 
defined with the correct Eq normalization factor required for the RG analysis [27]. 

The USSM particle content, in addition to the MSSM particles, includes a single 
extra scalar state, a new Abelian gauge boson and an additional neutral higgsino and 
gaugino state. The chargino sector remains unaltered, while the sfermion scalar potential 
receives additional D-terms. 



2.1 The Abelian gauge sector 

With two Abehan gauge factors, U(l)y and U(l)x, the two sectors can mix through the 
couphng of the kinetic parts [28] , 

>Cgauge = -^y" d^e [WyWy + WxWx + 2smxWYWx} , (2) 

where Wy and Wx are the corresponding chiral superfields. The gauge/gaugino part of 
the Lagrangian can be converted back to the canonical form by the GL(2,]R) transfor- 
mation from the original superfield basis Wy, Wx to the new one Wb, Wb' [28, 29]: 

Wy \ f 1 -tanx \ f Wb 



Wx J V l/cosx J yWB' J ' 
This transformation alters the U(l)y xU(l)x part of the covariant derivative to 

= d^ + zgYY,B^ + t{-gYYitsinX+ ^Qf)B'^ (4) 

cos X 

= d^ + ig,YiB^ + ig[Q[B'^, (5) 

where we introduced the notation gi — gY, g'l — gx/ cosx. We will also use g' — gi\/3/5 
for the low-energy (non-GUT normalized) hypercharge gauge coupling. 

With the above mixing matrix the hypercharge sector of the Standard Model is left 
unaltered, while the effective U(l)x charge is shifted from its original value Qf to 

Q[^Qf-^Y,t^nx. (6) 
9i 



5 



As a result of the kinetic mixing, new interactions among the gauge bosons and matter 
fields are generated even for matter fields with zero U(l)x charge. 

In the EgSSM the two U(l) gauge groups are automatically orthogonal at the GUT 
scale and the RG running effects give acceptable small mixing at the low scale [27, 30] 
providing in a natural way the general agreement between SM analyses and precision 
data [31]. Therefore in the reminder of the paper we will simply write Qi instead of Q[. 

After breaking the electroweak and U(l)x symmetries spontaneously due to non-zero 
vacuum expectation values of the iso-doublet and the iso-singlet Higgs fields, 

, sin/3 / \ , cos/3 / V \ 1 



^/2 \ v J' V2V0y' V2 

the Z, Z' mass matrix takes the form 

Mzz' -[ ^2 ^2 J , (8) 



where 



J2 _|_ „2 
T 



= ^iV^^' + ^2\ 2 (g^cos^/3-Q2sin^/3) 
Ml = g'y{Q\cos^P + Qlsin^P)+g'^Qlvl (9) 
We then diagonalise the mass matrix to give the mass eigenstates: 
Zi\ _ ( Z \ _ f cosazz' sinazz' \ f Z \ 

Z2 )~ \Z' )- \ -sinazz' cosazz' ) \ Z' ) ^'""^ 
where the resultant masses and ZZ' mixing angle are given by 

Ml,z, = ^ (^1 + ^1' T ^(M|-M|,f + 4A4^ (11) 

1 / 2A2 \ 

azz^ = -arctan(^ ^, (12) 

2.2 The Higgs sector 

In the charged sector it is convenient to introduce the G^, basis as: 

G- = cos 13 - H+* sin /3 (13) 
H+ = HJ* sin (3 + H+ cos (3 
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After the gauge symmetry breaking, two Goldstone modes from the original if„ and 
Hd doublets are eaten by fields leaving two physical charged Higgs bosons H^, with 
the mass 

where the trihnear coupling Ax is the soft-SUSY breaking counterpart of A, and the one- 
loop corrections A± are the same as in the MSSM [32] with the effective ji parameter 
given by 

In the CP-conserving model the CP-even and CP-odd scalar Higgs component fields 
do not mix. The CP-cvcn sector involves Reif°, Reif° and Re 5* fields. The 3x3 
mass matrix of the CP-even Higgs scalars Mjyg„ has been calculated to one-loop in 
Refs. [27, 33] in the field space basis /i, i7, S. This basis^ is rotated by an angle /3 with 
respect to the interaction basis. 
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\/2Re Hi \ = \ sin/? cos/? + wsin/? (15) 



The explicit form of Mg^g„ is given in Appendix A. It can be diagonalized by a 3 x 3 
orthogonal mixing matrix (O), i.e. 

J^2d^a, ^ O^Ml^^O (16) 

by going to the mass eigenstates basis 

{H,,H2,H,) = {h,H,N)0 (17) 

in which, by convention, mass eigenstates are ordered by mass, mij. < niHi+i- 
It will be convenient to introduce a mixing matrix O', 

(cos /? — sin /? \ 
sin/3 cos/? O, (18) 
ij 

that enters the Feynman rules. It is a superposition of two rotations in eqs. (15) and (17) 
and links the interaction eigenstates H^, H^, S directly to the CP-even mass eigenstates 
Hi, H2, H3. 



^Note that h, H are not the MSSM-like eigenstates. 
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The imaginary parts of the neutral components of the Higgs doublets and Higgs 
singlet compose the CP-odd sector of the model. In the field basis A, G, G' defined by 

v^Imi/° = G'cos/5 + (Acqs (/)-G" sine/)) sin/? 
V2Im//° = -G sin 13 + {A cos (t)-G'sm(t)) cos (3 
^ImS = Asin(/) + G"cos(/) (19) 

the massless pseudoscalar G, G' fields are absorbed to Z, Z' after the electroweak gauge 
symmetry breaking. The physical CP-odd Higgs boson A acquires mass 

rriA = . ^ , V IXea (20) 
sm20 

where tan0 = v sm.2(3 /2vs and the one-loop correction ^ea is given in Appendix A. 

Note that the Higgs sector of this model involves only one physical CP-odd pseu- 
doscalar as in the MSSM, since, unlike the NMSSM, the extra CP-odd state arising 
from the singlet is eaten by the Z' . However, there are three CP-even scalars, one more 
than in the MSSM, where the extra singlet state arises from the extra singlet as in 
the NMSSM. The characteristic Higgs mass spectrum in this model is governed by the 
value of A. For small values of A, say A < (71, the Higgs spectrum resembles that of the 
MSSM, with the heaviest CP-even Higgs scalar being predominantly composed of the 
singlet scalar state, and being approximately degenerate with the CP-odd pseudoscalar 
and the charged Higgs states when their masses exceed about 500 GeV. In this regime 
the hghtest CP-even Higgs scalar is Standard Model hke, and respects the MSSM mass 
bound. On the other hand, for large values of A, say A > 5^1, a viable Higgs mass 
spectrum only occurs for a very large CP-odd Higgs mass, say ruA ~ 2 — 3 TeV, with 
the heaviest CP-even Higgs scalar being non-singlet and degenerate with the the CP- 
odd and charged Higgs states. The second heaviest Higgs scalar is comprised mainly 
of the singlet state and is thus unobservable, while the lightest CP-even Higgs scalar is 
Standard Model like but may significantly exceed the MSSM bound. For more details 
concerning the Higgs sector see [27]. 



2.3 The neutralino sector 

The Lagrangian of the neutralino system follows from the superpotential in Eq. (1), com- 
plemented by the gaugino SU(2)£,, U(l)y and U(l)x mass terms of the soft-supersymmetry 
breaking electroweak Lagrangian: 

'^mi?"" = -\m2W''W'' -\myYY -]-MxXX - MyxYX ^h.c. (21) 

where the W"- (a = 1, 2, 3), Y and X are the (two-component) SU(2)l, U(l)y and U(l)x 
gaugino fields, and Mj {% — 2, X, y, YX') are the corresponding soft-SUSY breaking mass 
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parameters. After performing the transformation of gauge superfields to the gauge boson 
eigenstate basis, Eq.(3), the Lagrangian takes the form 

>^m^?"° = --M2W''W'' --MiBB --M[B'B' - MkBB' + h.c., (22) 



where 



M[ = 



Mx 2 sin y , ^ , ^ , 
cos^ X cos^ X 



-Myx + My tan^ x , Mk = 



YX 



cosx 



- My tan X, (23) 



and we introduce the conventional notation for the U(l) bino mass Mi = My. In 
parallel to the gauge kinetic mixing discussed in Sect. 2.1, the Abelian gaugino mixing 
mass parameter Myx is assumed small compared with the mass scales of the gaugino 
and higgsino fields. 

Notice that the gauge kinetic term mixing (and the corresponding soft-SUSY break- 
ing mass) can be a source of mass splitting between the B and B' gauginos in models 
with universal gaugino masses Mx — My. Since the mixing angle x must be small, as 
required by data [31], the splitting is very small. The splitting could be enhanced if 
additional U(l) gauge factors in the hidden sector were present that mix via the kinetic 
term with the visible sector.^ In our phenomenological analyses, therefore, we will con- 
sider two scenarios: (A) with M[ taken as a free parameter, independent from Mi; and 
(B) with M[ tight to Mi and M2 by a unification of gaugino masses at the GUT scale. 

After breaking the electroweak and U(l)x symmetries spontaneously the doublet 
higgsino mass /i and the doublet higgsino-singlet higgsino mixing /ix parameters are 
generated 



11 = X 



vs 
V2 



and iix = X 



x/2 



(24) 



The USSM neutral gaugino-higgsino mass matrix in a basis of two-component spinor 
fields ^ = {B, W^, H^, S, B'Y can be written in the following block matrix form 
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(25) 



Since the fields in the hidden sector are generally considered to be heavy enough and the hidden- 
visible mixing is expected to be small, their effect on the visible gauge sector can be negligible. Never- 
theless, the mass of the Abelian gaugino in the visible sector can obtain a substantial contribution, as 
advocated in [34]. 
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where the upper-left 4 x 4 is the neutral gaugino-higgsino mass matrix of the MSSM, the 
lower-right 2x2 corresponds to the new sector containing the singlet higgsino (singlino) 
and the new U(l)-gaugino B' that is orthogonal to the bino B, and off-diagonal 4x2 
describes the coupling of the two sectors via the ncutralino mass matrix {sp = sin (3, 
Cjj = cos/3, and sw,cw are the sine and cosine of the electroweak mixing angle 6*14^)- 
Notice the see-saw type structure of the new sector due to the absence of a diagonal 
mass parameter for the singlino ^S", which is in direct contrast to the NMSSM in which 
the cubic self-interaction generates a singlet mass term[10]. For the same reason, in the 
USSM the lightest neutralino can never be bino'-dominated. 

In general, the neutralino mass matrix M^o is a complex symmetric matrix. To 
transform this matrix to the diagonal form, we introduce a unitary 6x6 matrix N such 
that 

Xl = {B, Hd, Hu, S, B'), , (26) 

where the physical neutralino states x^, [A; = 1, 6] are ordered according to ascending 
absolute mass values. The eigenvalues of the above matrix can be of both signs; the 
negative signs are incorporated to the mixing matrix A^. Mathematically, this procedure 
of transforming a general complex symmetric matrix to the diagonal form with non- 
negative diagonal elements is called the Takagi diagonalization, or the singular value 
decomposition [18, 35]. Physically, the unitary matrix N determines the couplings of 
the mass-eigenstates xl to other particles. 

Although the complexity of neutralino sector increases dramatically by this extension 
as compared to the MSSM (which can be solved analytically), the structure remains 
transparent since, in fact, the original MSSM and the new degrees of freedom are coupled 
weakly. must be small by the requirement that the mixing of the U{l)x and U{1)y 
sectors satisfy experimental limits. The remaining off-diagonal terms are suppressed 
with respect to the corresponding block diagonal terms by a factor of v/vs- Since vs 
sets the mass of the Z', this results in vs being roughly an order of magnitude greater 
than V. Therefore in physically interesting case of weak couplings of both the MSSM 
higgsino doublets to the singlet higgsino and to the U(l)x gaugino, and the coupling 
of the U(l)y and U(l)x gaugino singlets, the remaining terms in the off-diagonal 4x2 
submatrix in Eq. (25) are small. Then, an approximate analytical solution can be found 
following a two-step diagonalization procedure given in Ref. [18]. In the first step the 
4x4 MSSM submatrix and the new 2x2 singlino-U(l)x gaugino submatrix At 2 
are separately diagonalised. In the second step a block-diagonalization removes the 
non-zero off-diagonal blocks while leaving the diagonal blocks approximately diagonal 
up to second order, due to the weak coupling of the two subsystems. 
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2.4 The sfermion sector 

As explained in the Introduction, we assume the exotic squarks to be substantially 
heavier than the MSSM fields. However the structure of the MSSM squarks gets modified 
by the presence of extra U(l)x- Both the squarks and sleptons are important to our 
analysis and so we briefly describe the new ingredients in the sfermion mass matrix 
(neglecting the possibility of flavor and CP violation) 

f \ m/(A/ -//(tan/3) ^^f) m^^ + mj + Aj, I ^ ' 

where rap, nij are the sfermion soft-supersymmetry breaking parameters for the quark 
and lepton doublets F — Q,L and singlets / = U^,D'^,E'- and ^/ is the trilinear 
coupling, while m/ is the corresponding fermion mass and the D-terms receive additional 
U(l)x terms 

A_^- = M| cos 2(3{lj - e/s^) + lg[^Qf [v^ {Qi cos^ p + Q^ sin^ (3) + Qgvl] (28) 

where Ij and e/ are the weak isospin and electric charge and the U(l)x charges Qf are 
for the left fields. Explicitly, we have for squarks 

= M|cos2/3(|-|s^) + ^y;2[-T;2(2sin2/3 + 3cos2/3) + 5T;|] 

As* = M| cos 2/3 + ^g[^[-v^{2 sin^ /3 + 3 cos^ /3) + Bv^] 

Aj = M|cos2/3(-i + |s^) + g^^;2[_^2^2sin2/3 + 3cos2/3) + 5^;|] 

Aj, = M|cos2/3(-is2^) + ^^;2[-v2(2sin2/3 + 3cos'/3) + 54] (29) 

and for sleptons 

Ap = Mlcos2(3{^) + ^g[^[-v\2sm^(3 + 3cos^(3) + 5vl] 
Ae = Mlcos2f3{-l + sl,) + ^g[^[-v\2sm^(3 + 3cos'(3) + 5vl] 
Ae* = M| cos 2/3 (-s^) + ^^;2[-v2(2 sin' /3 + 3 cos2/3) + 5v|] (30) 

Note that here g[ is the GUT normahzed U{l)x gauge coupling analogous to the GUT 
normalized hypercharge gauge coupling gi in the MSSM. 

The diagonal form of the sfermion mass matrix is obtained, as usual, by a 2x2 
rotation in the LR plane 

^2diag ^ UjMjUf (31) 

and the mass eigenstates are defined according to 

/l ^ _ 7-rt f fL 



f2 J '\fR 

with the convention that m; < m; . 

h — h 
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3 Calculating the relic density 



The calculation of the neutralino LSP relic density in the MSSM is well known [1, 2] and 
has been widely studied in the general MSSM [36] and the constrained MSSM [37]. The 
calculation of the relic density in the NMSSM has also been extensively studied [38]. 
The differences between the MSSM relic density calculation and the USSM calculation 
arise through the extension of the particle spectrum and through the new interactions 
that are introduced. We have implemented all new interactions into the micrOMEGAs[25] 
code using LanHep[24] to generate the feynman rules. MicrOMEGAs takes full account of 
all annihilation and coannihilation processes and calculates their effect whenever they 
are relevant. Nevetheless, from the form of these alterations we would hke to make some 
general observations before we go on to consider the details of the calculations. 

The USSM extends the neutralino sector by adding two new states to the spectrum: 
the bino' and singlino components. This results in two extra neutralinos. However 
for the relic density calculation we are only interested in the lightest neutralinos, so 
the primary effect will be through the magnitude of the singlino and bino' components 
in the lightest neutralino. In what follows we will be interested in the scenarios in 
which the lightest neutralino has a significant singlino component and a small but non- 
zero bino' component. Therefore it is informative to consider the general form of the 
interactions that arise from the singlino and bino' components of the lightest neutralino 
before considering specific diagrams. 

The bino' component is always sub dominant to the singlino component due to the 
see-saw structure of the extra 2x2 S / B' sector of the neutralino mass matrix in Eq. 25. 
The form of the interactions that arise from the inclusion of the bino' component closely 
mirror those of the bino component, except for the different coupling constant and 
charges under the new U{l)x- 

The singlino component is another matter. It gives rise to a new type of neutralino 
interaction from the \S HuHd term in the superpotential that will be seen to dominate 
the annihilation processes of neutralinos with significant singlino components. This term 
means that if the lightest neutralino has significant singlino and higgsino components 
then it will couple strongly to Higgs bosons with a significant or component, 
usually the lighter Higgs bosons, Hi -2 and A in the spectrum. Moreover, the absence 
of the singlet cubic term S^, in contrast the the NMSSM, implies that the singlino- 
dominated LSP needs an admixture of MSSM higgsinos to annihilate to Higgs bosons. 

On the other hand, the singlino component does not interact with the SU{2) or 
U{1)y gauginos. Therefore a significant singlino component in the lightest neutralino 
will suppress annihilations to W or Zi bosons. 

Finally, there is no coupling of the singlino component to fermions. Thus a significant 
singlino component in the lightest neutralino will also suppress annihilation to fermions. 

Having noted these general features we will now consider the specific behavior of the 
different annihilation diagrams. 
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3.1 t-channel diagrams 

• Gauge boson final states 




Figure 1: The t-channel annihilation processes for a neutralino to final states containing 
gauge bosons. 

Fig. 1 shows the t-channel diagrams available for annihilation of neutralinos to gauge 
bosons. The XiX'jZi vertex is given in Eq. (B.l). Note that the coupling of neutralinos 
to the Z component of the Zi state is precisely that of the MSSM XiX^^ coupling. As 
the Z component dominates the Zi state, a singlino dominated LSP will not annihilate 
strongly to Zi bosons. 

In contrast there is a strong coupling from the MSSM-higgsino components as well 
as the singlino component to the Z' component of the Z^ state. Notice also that the 
MSSM-higgsino components of the LSP enter with the same sign in the coupling to 
the Z', unlike in the coupling to the Z, where they tend to cancel each other. As the 
Z2 boson is dominantly Z' any LSP with a non-zero higgsino or singlino fraction will 
annihilate to Z2 bosons when such a final state is kinematically allowed. Unfortunately 
the Z2 is required to be heavy by experimental limits, so annihilation of the lightest 
neutralinos to final states involving one Z2 is hard to achieve and annihilation to two 
Z2 bosons is impossible. 

The second diagram of Fig. 1 shows the t-channel annihilation to final states. 
Eqs. (B.2) and (B.3) give the relevant coupling and show that the singlino and bino' 
components do not couple to the wino component of charginos or to the bosons. 
This means that a large singlino or bino' component in the LSP will suppress annihilation 
to bosons in the final state. 

• Higgs boson final states 

Fig. 2 shows the available t-channel processes for the annihilation of neutralinos 
to final state Higgs bosons. Due to the XS HuHd term in the superpotential and the 
D-terms there are significant differences between these diagrams in the USSM and the 
MSSM. The XiXiHj vertex given in Eq. (B.6) is the relevant vertex in this first diagram. 
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Figure 2: The t-channel annihilation processes for a neutralino to final states involving 
scalar Higgs bosons, pseudoscalar Higgs bosons or charged Higgs bosons respectively. 

First note that the bino' component of one neutrahno couples with the higgsino 
component of the other and the component of the final state Higgs boson in the 
same way as the equivalent coupling of the bino or wino components. In addition there 
is an extra term which couples the bino' component of one neutralino to the singlino 
component of the other and to the singlet component of the Higgs boson in the final 
state. This means that if the lightest neutralino is dominantly singlino, it will annihilate 
to final state Higgs bosons with a significant singlet component through the exchange of 
a neutralino with a significant bino' component in the t-channel. Unfortunately these 
processes are disfavored for the same reason as annihilation to final states containing a 
Z2. The Higgs boson with a significant singlet component will have a mass comparable 
to the Z2 boson and thus a final state with two such Higgs bosons will be impossible 
and even one will often be kinematically ruled out. 

Of more interest is the term in this vertex that couples a singlino component of one 
neutralino to a higgsino component of the other neutralino and the iJ^^^ components of 
the Higgs boson with a strength A. If the lightest neutralino is dominantly singlino then 
two LSPs can exchange a dominantly higgsino neutralino in the t-channel to produce 
two Higgs bosons in the final state. This is a channel that is always present if the 
lightest neutralinos are heavy enough to produce two light Higgs bosons in the final 
state. Obviously, if both Hi and H2 are lighter than the hghtest neutrahno then there 
will be more available channels. As the singlino couples predominantly to Higgs states, 
this channel provides the strongest annihilation mechanism for a neutralino with a large 
singlino component. This amplitude will be maximised for three degenerate mixed 
state neutralinos with strong higgsino and singlino components that are heavier in mass 
than the hghtest two Higgs states. The addition of this vertex also allows for a new 
annihilation process for a dominantly higgsino neutralino through the exchange of a 
t-channel neutralino with a substantial singlino component. 

The middle diagram of Fig. 2 shows the annihilation to final state pseudoscalar 
Higgs bosons. The relevant vertex is given in Eq. (B.IO). The first line gives the 
familiar MSSM vertex for the coupling of a S or I^ component of a neutralino to a 
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higgsino component and a pseudoscalar Higgs. This is modified by an overall factor of 
cos (j) which determines the magnitude of the MSSM-like components of the pseudoscalar 
Higgs over the singlet contribution. As sin0 ^ v cos (3/vs, the suppression from cos(j) 
is small. This is the same as saying that the pseudoscalar Higgs generally only has a 
very small singlet component. The analogue of the W, B interaction terms appears for 
the bino'. The bino' component also couples to the singhno component of the second 
neutralino and the singlet component of the final state pseudoscalar Higgs. This term is 
sin<;/) suppressed due to the small singlet component of the A bosons in the final state. 
These interactions determine the strength of the annihilation of a dominantly gaugino 
LSP to pseudoscalar Higgs bosons through the exchange of a dominantly higgsino (or 
singlino) neutralino. 

More interesting contributions come from the \S HuHd term in the superpotential. 
These provide a AHuHd coupling, albeit suppressed by a factor of sin0. Such a coupling 
does not appear in the MSSM. There is also a term that couples ASH^^d with no sin 
suppression. Once again this produces a strong annihilation channel for a neutralino 
with a substantial singlino component through t-channel neutralino exchange where the 
neutralino exchanged in the t-channel must have a significant higgsino component. This 
is the analogue of the process we discussed in some detail for the scalar Higgs final 
states and will, kinematics allowing, give a strong annihilation channel for a dominantly 
singlino neutralino as long as there is a light neutralino in the spectrum with a substantial 
higgsino component to be exchanged in the t-channel. 

The final diagram of Fig. 2 shows annihilation to charged Higgs boson final states. 
The relevant vertex is given in Eq. (B.14). The vertex includes a B' interaction that 
parallels the familiar B and W interactions to the higgsino component of the chargino 
and a charged Higgs boson. There is also a term that arises from the H^Hd super- 
potential term. This allows for a neutralino with a substantial singlino component to 
annihilate to charged Higgs bosons via t-channel chargino exchange as long as there are 
light charginos with a significant higgsino component and the final state charged Higgs 
bosons are kinematically allowed. In contrast to the previous two diagrams, this one 
does not add an extra annihilation channel for a dominantly higgsino neutralino. In the 
first two diagrams there is the new possibility in which a dominantly singlino neutralino 
is exchanged in the t-channel. In the third diagram there is no such process as there is 
no singlino component in the charginos. 

From an analysis of the processes with Higgs bosons in the final state we see that 
there will be a strong annihilation cross-section for a neutralino with a large singlino 
component to light Higgs bosons if there is a light neutralino with a substantial higgsino 
component in the spectrum and the Higgs boson final states are kinematically allowed. 
We also note that the \S HuHd allows for new couplings between neutralinos and Higgs 
bosons that will alter the annihilation of dominantly higgsino neutralinos with respect 
to their behavior in the MSSM. 
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Figure 3: The t- channel annihilation process for neutralinos to fermions. 

• Mixed boson final states 

It is quite possible to have an unmatched pair of bosons in the final state of a t- 
channel annihilation diagram. We do not need to go through the details of all possible 
diagrams. Instead we just note that a neutralino with a significant singlino component 
will dominantly annihilate to final states made up of Higgs bosons. The strength of such 
channels will depend upon the size of the singlino component in the lightest neutralino, 
the mass of the neutralinos with substantial higgsino components that will be exchanged 
in the t-channel, and the mass of the final state Higgs bosons. 

• Fermion final states 

Finally we consider the t-channel annihilation diagram to final state fermions through 
the diagram given in Fig. 3. The squark vertices are given in Eqs. (B.18) and (B.19). 
The couplings of the bino and wino components of the neutralino are the same as in the 
MSSM. Note that there is an extra couphng of the bino' component of the neutrahno to 
the squark-quark pair that is of the same order of magnitude as for the B. As the bino' 
is only ever a subdominant component of the lightest neutralino, and as the annihilation 
to fermions is relatively weak in the first place, we can expect that interactions of this 
form will have little impact on the annihilation cross-section. However, if the lightest 
neutrahno is too light to annihilate to final state Higgs bosons, this channel will remain 
open and can dominate though it will give a relic density well in excess of that measured 
by WMAP. 

3.2 s- channel diagrams 

Fig. 4 shows the possible s-channel processes available for the annihilation of a pair 
of neutralinos. The first diagram shows the annihilation through and intermediate Zi 
gauge boson. The relevant coupling of two neutralinos to a Zi is given in Eq. (B.l). 
As before we note that the singlino component of the neutralino only couples to the 
Z' component of the Z^ gauge boson. This means that if the lightest neutrahno has 
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Figure 4: The annihilation processes for a neutralino through s- channel Higgs and Zi 
bosons, where we do not specify the precise particles in the final state. 

a significant singlino component, then annihilations through an s-channel Zi will be 
suppressed as the Z' component of the Zi is required to be very small. On the other 
hand, the Z2 has a large Z' component. Therefore a lightest neutralino with a substantial 
singlino component will annihilate through an s-channel Z2. 

The second diagram shows the annihilation of ncutralinos through an s-channel 
scalar or pseudoscalar Higgs boson. The relevant couplings are given in Eq. (B.6) and 
Eq. (B.IO) respectively. Consider first the case in which the s-channel Higgs boson is 
dominantly composed of the singlet Higgs. Since the bino or wino components of the 
hghtest neutralino only couple to the non-singlet Higgs components of the s-channel 
Higgs a light neutralino that is dominantly wino or bino will not annihilate strongly 
through a dominantly singlet Higgs. However there is a coupling of the singlet com- 
ponent of the Higgs boson to the higgsino components of the lightest neutralino. This 
provides a strong channel when on-resonance for annihilation of a light neutralino with 
a large higgsino component. There is also a strong couphng if the lightest neutralino 
has significant bino' and singlino components. Thus wc expect a light neutralino with 
strong mixing between higgsino, singlino and bino' terms to annihilate strongly through 
s-channel heavy Higgs exchange where the heavy higgs has a large singlet Higgs compo- 
nent. 

If the s-channel Higgs boson does not have a large singlet Higgs component then the 
story is somewhat different. In this case the light neutralino needs to have a significant 
higgsino fraction along with a substantial contribution from one of the other non-higgsino 
states. This situation is mirrored in the case of the pseudoscalar Higgs. 

Prom this we see that we have a new annihilation channel for neutralinos with a 
significant higgsino fraction through a dominantly singlet Higgs in the s-channel. We also 
see that a light neutralino with a substantial singlino-higgsino mixture will annihilate 
strongly through the whole range of s-channel Higgs exchange processes. 
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3.3 Coannihilation 

As well as the annihilation of two identical neutralinos, it is often the case that coannhi- 
lation between the LSP and the NLSP (and sometimes even heavier states) can be 
important. This process is normally important for a dominantly MSSM-higgsino or 
wino neutralino LSP. In these situations there is an automatic near degeneracy in the 
mass of the lightest neutralino with the mass of the lightest chargino and, in the case 
of the higgsinos, also with the next-to-lightest neutralino. A dominantly singlino LSP 
does not have an automatic degeneracy with other states. However, it is possible for a 
singhno neutralino to be exactly degenerate with other states - something that does not 
happen in the MSSM due to the signs of the terms in the neutralino mixing matrix. In 
these cases wc would expect the effect of coannihilation to be important. 

Therefore we expect coannihilation processes to only be significant in regions of the 
parameter space where we move from one type of LSP to another as this indicates a 
degeneracy in the mass of the LSP and NLSP. We also expect to see the standard large 
coannihilation contributions for a predominantly MSSM-like higgsino LSP or predomi- 
nantly wino LSP. 

4 Elastic scattering of neutralinos from nuclei 

The direct cold dark matter search experiments, such as DAMA/LIBRA, CDMS, ZEPLIN, 
EDELWEISS, CRESST, XENON, WARP [39], aim at detecting dark matter particles 
through their elastic scattering with nuclei. This is complementary to indirect detection 
efforts, such as GLAST, EGRET, H.E.S.S. [40], which attempt to observe the annihila- 
tion products of dark matter particles trapped in celestial bodies. 

Since we assume the LSP to be the hghtest neutrahno x^, we consider the elastic 
scattering of the hghtest neutrahno from nuclei. The elastic scattering is mediated by 
the t-channel Zi and Higgs exchange, as well as the s-channcl squark qj exchange, 
as depicted in Fig. 5 for x5 'i scattering. There are also important contributions from 
interactions of neutralinos with gluons at one loop [42, 43]. 



Figure 5: Diagrams contributing to the lightest neutralino scattering from a quark. 
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The extended particle content and new couplings present in the USSM model have 
also a direct effect on the elastic cross section calculations, as discussed in the previous 

chapter. 

The clastic cross section for ncutralino scattering from a nucleus can be broken into 
a spin-independent (SI) and a spin-dependent (SD) part, 

(7 = CTgi (7gi3 , (33) 

each of which can be expressed in terms of the elastic scattering of neutralino from 
individual nucleons in the nuclei. In the limit of zero-momentum transfer they can be 
written as [41] 

A 2 

a^^ = ^[ZU + {A-Z)U]\ (34) 

TT 

3977?^ 

asD = ^G|J(J+l) A^ (35) 

71 

where Z and A are atomic number and mass of the nucleus, J is the total nucleus 
angular momentum and is the reduced neutralino-nucleus mass. Note that the spin- 
independent part benefits from coherent effect of the scalar couplings, which leads to 
cross section and rates proportional to the square of the atomic mass of the target nuclei. 
The spin-dependent quantity A is given by 



, 1 



g=u,d,s ^ q=u,a,s 



(36) 



where (Sp) and (S'„) are the expectation values of the spin content of the proton and 
neutron group in the nucleus, while and are the quark spin content of the proton 
and neutron, respectively. 

For the spin-independent part, the effective couplings of the LSP neutrahno to proton 
and neutron fp and are more complicated. In the limit of m^o <^ rriq and ruq <^ rriq, 
which we will later confine to, they simplify and can be approximated as: 

mp,n ^ ruq 27 ■'-'G ^ m, ^ ^ 

q=u,d,s ^ q=c,o,t ^ 

The first term in Eq.(37) corresponds to interactions with the quarks in the target nuclei, 
while the second term corresponds to interactions with the gluons in the target through 
a quark/squark loop diagram, and 

q=u,d,s 
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Finally, the effective Lagrangian for elastic scattering of neutralinos from quarks in the 
non-rclativistic limit can be written as a sum of axial- vector (spin-dependent) and scalar 
(spin-independent) terms: 

£eff = Ag ixiYlbXi) (Qlnlbq) + Bg (xiXi) (qq) (39) 
The effective couplings Ag and Bg are given by: 
r,2 I R'-^-P -I- I R«-R|^ n 

[gi|iVi3|' + g2|iVi4|' + g.|iVi5i'] (Qq + q^) m 



4m|, 



=1,2 



hg ^ Re(Gfe) + Re{G'^) + RejGl) f 0[, foTq = d,s,b 
2V5^ \ Oafe ior q = u,c,t 

In this expressions we have neglected a small Z-Z' mixing. 

The first terms in both effective couplings come from squark exchange diagrams. 
The neutralino-squark-quark couplings 5*^, 5*^ are given in Appendix B. As seen in 
Eqs. (B.20,B.21), they receive a contribution form the bino' component Niq. 

The second and the third terms in (40) come from the Z and Z' exchanges, re- 
spectively, where the latter contains a term due to the singlino component, A^is. The 
second term in the form factor Bg receives contributions from three scalar Higgs boson 
exchanges. Each contains an MSSM-like term, Gk, as well as the new terms G'^ and G'^, 
(A; = 1,2, 3) 

Gk = 92iNu-twNn){N^^O'^,-NnO[k) 

G', = -2g[N,e{QiN,sO[, + Q2NuO'^,+QsN,,0',,) 

G'l = ^\[N,,{N,^0'^^ + N^^O[^) + N,^N,,0'^^] (42) 

The G'^. piece is generated by the g[B'{HiHi + SS) couplings from the extra U(l)x 
D-terms, while the G'^ is induced by the XHi{SHj + HjS) couplings (here we follow the 
conventions and notations of Ref. [44], properly extended to the USSM model [45]). 

5 Results 

Now that we have introduced the model we move on to study the details of the dark 
matter phenomenology within the USSM parameter space. 
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5.1 Defining a parameter range 

Before we study the phenomenology we need to define the parameter range we are 
interested in. The USSM extends the number of free parameters over those in the 
MSSM by the set: 

M(, ^i, A, Ax, vs. 

These parameters are constrained by a number of different factors. 

We fix g'^ = Qi as we wish to maintain gauge couphng unification and the two U{1) 
gauge couplings run with identical RGEs. 

The parameters vs, A, Ax appear in the determination of particle masses. Therefore 
we determine these by setting the corresponding masses. First of all, we wish to keep vs 
low to maximise the region of parameter space in which there is a light singlino/bino' 
LSP. If M[ = then there are two degenerate singlino/bino' neutralinos with a mass 
Q'sQ'iVs- However, we do not have the freedom to set Vs arbitrarily low since from 
Eqs. (9) we see that low vs would require a light Z2 mass and a large Z-Z' mixing 
incompatible with the LEP and Tevatron hmits. Adopting 

rris = g[vs = 1200 GeV, (43) 

together with assumed tan/3 = 5, gives = 949 GcV and sinzz' = 3 10~^ which is 
consistent with current constraints. We use this to set the magnitude of the vs in all 
that follows. 

With Vs set, our choice of A will set the size of // through the relation 

,= A^. (44) 

Note that A is a coupling and so cannot be too large. An upper limit on A < 0.7 at 
a given value of vs results in a corresponding maximum value on /i, and consequently 
/J, < ms,z'- As a result, m^o < 777,^/^5 will always be satisfied which has important 
implications for the available dark matter annihilation channels. It also justifies our 
earlier claim that there will always be light charginos and higgsinos in the spectrum if 
the Z' mass is low. 

We set the size of Ax by setting the mass of the pseudoscalar Higgs. Prom Eq. (20) 
we see that once the VeV of S and A have been set, the mass of only depends upon 
tan/3 and Ax- As we are keeping tan/3 fixed, we can use Ax to set the psuedoscalar 
Higgs mass. 

The familiar MSSM parameters are also relevant to the details of both the relic 
density calculation and the direct detection phenomenology. The most important pa- 
rameters are those that appear in the neutralino mass matrix - Mi and M2. We keep the 
ratio Ml : M2 = 1 : 2 for simplicity, but there are as many ways to break this relation 
in the USSM as in the MSSM. 
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Finally we must set M[. In what follows we take M[ as a free variable and scan 
over a range of values. In our first study we take M[ to be independent of the other 
gaugino masses, as in the study of Rcf. [18] where the collider phenomenology has been 
discussed. This will complement Ref. [18] with the dark matter calculations. On the 
other hand, it is also interesting to consider a scenario in which soft SUSY breaking 
gaugino masses are unified, namely Mi : M[ : M2 = 1:1:2 and we do this in our 
second scenario. This allows us to organize our studies in the following way: 

• scenario A: M[ arbitrary ; 

• scenario B: unified gaugino masses M[ = Mi = M2/2 . 

To calculate the relic density we need to set the rest of the particle spectrum. To do 
this we fix the pseudoscalar Higgs mass = 500 GcV and for sfermion masses we take 
a common mass of mQ^u,d,L,e = 800 GeV, and a common trilinear coupling A = 1 TeV, 
while the gluino mass is determined assuming unified gaugino masses at the GUT scale. 
We have set the squarks and sleptons to be heavy as this allows for a clearer analysis of 
the annihilation properties of the neutralinos. 

For the direct dark matter searches, there are large uncertainties in the spin-dependent 
and spin-independent elastic cross section calculations due to the poor knowledge of the 
quark spin content of the nucleon and quark masses and hadronic matrix elements. 
These uncertainties have recently been discussed in Ref. [46], from where we calculate 
the central values of f^'^: 

fl = 0.027, fl = 0.039, =0.36 

J^d '-a 

r = 0.0216, r = 0.049, r = 0.36 (45) 



and AP'": 



= +0.84, A^ = -0.43, A^ = -0.09 
a;] = -0.43, A2 = +0.84, A^ = -0.09 (46) 



5.2 Scenario A: M[ arbitrary 

In this scenario we take M[ as an arbitrary parameter with the MSSM gaugino param- 
eters fixed at // = 300 GeV, Mi = M2/2 = 750 GeV. 



5.2.1 Mass spectrum 

With these parameters we calculate the resulting mass spectrum at a given value of M(. 
The mass spectrum for M[ = GeV is shown in Fig. 6. In the Higgs sector we have 
a light Higgs at 127 GeV, a heavier scalar, pseudoscalar and charged Higgses around 
500 GeV and a dominantly singlet Higgs at 949 GeV. Sfermions are located between 750 
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Parameter 


Value 


Ml 


750 GeV 


M2 


1500 GeV 


/i 


300 GeV 


M[ 


0-20 TeV 


(S) 


2607.61 GeV 


A 


0.163 


Ax 


160 GeV 



Table 2: The parameters taken for the neutralino sector in the scan with /i = 300 GeV, 
ruA = 500 GeV, tan /3 = 5. 
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Figure 6: T/ie mass spectrum for M[ = GeV 

to 950 GeV. The chargino sector consists of a higgsino-like chargino around 300 GeV 
and a wino-like chargino at 1500 GeV. Since the mixing between the MSSM-like and 
the bino'/singlino at M[ = GeV is numerically small, the spectrum of neutralinos can 
qualitatively be understood by separately diagonalizing the 4x4 and 2x2 neutralino mass 
sub-matrices. Thus to a good approximation we have (according to ascending (absolute) 
masses for M[ = GeV) a pair of nearly degenerate, maximally mixed MSSM higgsinos 
at 300 GeV (first two states), an MSSM bino at 750 GeV (the third), a pair of nearly 
degenerate, maximally mixed singlino/bino' neutralinos at 949 GeV (the fourth and the 
fifth) and an MSSM wino at 1500 GeV (the sixth state). 

To understand the change of neutralino masses and of their composition as a function 
of M[ it is instructive the follow their analytic evolution as M( is turned on. The see-saw 
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Figure 7: The neutralino mass spectrum for varying M[, for the parameter choices in 
table. 2. The right panel is a magnified part of the left one. 

structure of the 2x2 singlino/bino' submatrix forces the two nearly degenerate, mixed 
singhno/bino' states to move apart: the hghter one (the fourth) gets hghter, and the 
mass of the other (the fifth one) heavier as M[ increases. The MSSM-hke states do 
not evolve much, unless the mass of one of the new states comes close to one of the 
MSSM, where a strong mixing may occur. For the mixing to be important not only 
the (absolute) masses must come close, but also the mass-eigenstates must belong to 
eigenvalues of the same sign. It is obvious from the see-saw structure that the heavier 
singlino/bino' state (the fifth one) belongs to the positive and the lighter (the fourth) to 
the one negative eigenvalue. Similarly the lighter of the two nearly degenerate MSSM 
higgsinos (the first state) belongs to the positive, and the other (the second) to negative 
eigenvalue. 

As M[ increases the (absolute) mass of the fourth state gets closer to the third, 
however they do not mix since they belong to eigenvalues of opposite sign. In left panel 
of Fig. 7 the lines representing these two states pass each other at M[ ~ 450 GeV. The 
bino, which is the third state according to the mass ordering below 450 GeV, becomes 
the fourth one when M[ passes 450 GeV. On the other hand when M[ approaches 900 
GeV and the mass of the fifth state gets close to the sixth one, strong mixing occurs 
between these states - the two lines representing these states in Fig. 7 "repel" each 
other. The heaviest neutralino smoothly changes its character from the MSSM wino 
to the singlino/bino' when M[ passes the cross-over zone near 900 GeV. Even more 
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interesting feature occurs when M[ approaches 2500 GeV, as illustrated in the right 
panel of Fig. 7 - a magnified part of the left panel. The singlino/bino' state belonging 
to the negative eigenvalue (which is now the third state according to mass ordering) 
mixes strongly with the second one. It does not mix with the first one since these 
states belong to eigenvalues of opposite sign. As a result of the mixing the mass of the 
second state is pushed down and below the lightest one for M[ above ~ 2.6 TeV. Thus 
the LSP discontinuously changes its character from being mainly higgsino to mainly 
singlino/bino' when M[ passes the cross-over zone near 2.6 TeV. For higher M[ values 
the LSP becomes dominantly singlino. This behavior will be important to understand 
discontinuities in plots to follow. 

5.2.2 Relic density 
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Figure 8: The relic density across varying M[, for n = 300 GeV, tra = 500 GeV and 
tan/5 = 5. The red lines show the 2(T measurement of the QcoMh"^ by WMAP-5. The 
green line shows the approximate MSSM higgsino relic density for ^ = 300 GeV. 

Having set the masses, we vary M( and calculate the relic density. The resulting 
values for the relic density are plotted in Fig. 8. Before dealing with the specific channels 
that give rise to the different features, we make some general points. Firstly, as /i = 
300 GeV and mA,H,H± ~ 500 GeV it is never possible for a pair of neutralinos to 
annihilate to a pair of pseudoscalar Higgs bosons, medium mass Higgs bosons or charged 
Higgs bosons in the final state. Secondly, as the squarks and sleptons are significantly 
more massive than the mass of the LSP, they do not contribute significantly to the 
annihilation cross-section except where noted below. 
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In the range < M( < 2.5 TeV the LSP is predominantly composed of MSSM- 
higgsino and gives a reUc density of the same order of magnitude as an MSSM-higgsino. 
At Ml = 2.57 TeV the LSP becomes dominantly singhno, as shown by the cross-over of 
the mass hncs in Fig. 7. As M[ increases, the singhno component of the LSP increases 
steadily. This decreases the strength of the Xi — Xi coannihilation that dominates the 
annihilation amplitude for a predominantly MSSM-higgsino LSP. As a result we might 
expect the value of flcDMh'^ to increase noticably before M[ = 2.57 TeV. However, as 
M[ approaches 2.57 TeV the mass splitting between Xi and X2 decreases. This increases 
the amplitude for Xi ~ X2 coannihilation. This increase compensates the drop in the 
neutralino-chargino coannihilation and results in an almost fiat value of flcDMh"^ up to 
M[ = 2.57 TeV. 

Above M[ — 2.57 TeV the mass splitting between the lightest neutralino and the 
lightest chargino and next to lightest neutralino increases steadily. This quickly turns off 
any coannihilation processes. At the same time, the singlino component of the lightest 
neutralino increases quickly. This steadily reduces the amplitude of x5~Xi annihilations. 
As a result of the combination of these two effects there is a sharp rise in the relic density 
above M[ = 2.57 TeV. 

At M[ ^ 3 TeV we see a sharp dip in the value of QcDMh'^ caused by the pscudoscalar 
Higgs s-channel resonance. Just below M[ = 5 TeV we see a sharp jump in the relic 
density as the LSP drops below the top mass, ruling out processes of the form x^Xi ~^ 
H* ti. By M{ = 5 TeV the LSP is 94% singlino with a 3% bino' admixture and a 
2% higgsino admixture. This, combined with the mass sphtting between the higgsinos 
and the singlino LSP, suppresses the annihilation of the singlino resulting in a relic 
density well above the measured value. At this point the dominant annihilation channel 
is to b, b through off-shell s-channel Higgs production, with a subdominant contribution 
from t-channel higgsino exchange to final state light Higgs bosons. A small kink in the 
relic density profile at M[ 11 TeV is the point at which the singlino becomes lighter 
than the light Higgs boson and final states with two Higgs bosons become kinematically 
disallowed. The dip at M[ = 14 TeV is the light Higgs resonance and the dip at 
M[ = 20 TeV is the Zi resonance. 

Here we have seen that the dominant annihilation channels of the singlino - through 
t-channel higgsino exchange and through s-channel Higgs production - are not strong 
enough to give a relic density in agreement with the measured value. The exception 
is when the singlino is mixed with a higgsino state. This enhances the annihilation 
through s-channel Higgs production as the neutralino-neutralino- Higgs vertices require a 
non-zero higgsino contribution. It also enhances annihilation through t-channel higgsino 
exchange as the higgsinos are fighter. 

The fact that we find a large relic density for a singlino LSP is partly down to our 
choice of parameters. Singlino dark matter dominantly annihilates to Higgs bosons, and 
with the parameters chosen above all but the lightest Higgs boson arc excluded from 
the final state by kinematics and s-channel processes are similarly suppressed by the 
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large masses. This would not be the case if we were to take mA,H,H± < We can do 
this by either lowering Ax or increasing A. Raising A also has the effect of increasing 
the coupling strength of the relevant vertices for singlino annihilation. We will discuss 
these effects further in scenario B. 



5.2.3 Direct detection 




Figure 9: The elastic spin-independent (left) and spin- dependent (right) LSP-proton 
cross section as a function of M[ in scenario A. 

Let us now turn to the direct DM detection analysis. In Fig. 9 the spin-independent 
as well as spin-dependent elastic cross section of the lightest neutralino on a proton is 
shown as a function of M[. We restrict the range of M[ to 2-5 TeV, since beyond this 
range the cross section changes monotonically. 

To understand the M[ behavior, we refer to Fig. 7. Up to M[ ~ 2.5 TeV the lightest 
neutralino is almost a pure MSSM higgsino. As a result its couplings do not depend 
on M[ and the scattering cross sections are practically determined by the MSSM-like 
terms Gk- Both the SI and SD cross sections are almost equal to the MSSM result with 
corresponding parameters. 

The discontinuity in the cross sections around 2.5 TeV is related to the sudden 
change of the nature of the LSP. As the M[ parameter increases, the mixing between 
the third and the second states pushes the latter below the lightest one (right panel 
of Fig. 7). The nature of the LSP therefore changes discontinuously from one of the 
MSSM-like higgsinos to the other higgsino state which at the same time acquires an 
increasing singlino component. 

The reduction of the spin-independent cross section (left panel) can be understood 
by reahzing that the elastic cross section of the second-lightest state (according to mass 
ordering below M[ = 2.5 TeV) on the proton is more than an order of magnitude smaller 
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than that for the hghtest one. When it becomes the LSP (for M[ > 2.5 TeV) the SI 
cross section drops significantly. As the singhno component of the LSP increases with 
M[ the G'l. factors, which are sensitive to both the singhno and the higgsino components 
- viz. Eq. (42), become responsible for the rise of the cross section. With further increase 
of M( the LSP becomes almost a pure singhno which explains a steady fall of the cross 
section. 

The spin-dependent cross section is dominated by the gauge boson exchange diagram. 
The Z coupling to the lightest neutralino is controlled by the combination C34 = |iVi3p — 
|iVi4p of neutralino mixing matrix elements. For low M[ the lightest neutralino is 
almost a perfect mixture of and for which these elements almost entirely cancel 
resulting in a small value of C34. As M[ increases the singhno forces the second-lightest 
state to become the lightest (flipping the sign of the coupling) and upsets this delicate 
cancelation. As a result, the cross section increases by a factor 6 and then starts to fall 
as the LSP becomes dominantly a pure singhno state. 

5.3 Scenario B: Mi = M[ = M2/2 

In the previous subsection we have considered the phenomenology of the USSM with 
non-universal Mi and M[. In this section we will consider the scenario in which gaugino 
masses are unified at the GUT scale implying the ratio Mi : M[ : M2 = 1 : 1 : 2 at 
the electroweak scale. We will vary M[ (together with other gaugino masses) as before 
and consider the behavior of both the relic density and the direct detection behavior. 
Motivated by the remarks at the end of Subsection 5.2.2 we also increase the value of /i 
parameter by a factor of 2, i.e. we take = 600 GeV. This is achieved by doubling the 
size of A. 

5.3.1 Mass spectrum 

Again to understand qualitatively the neutralino mixing pattern we start the discus- 
sion with M[ — 0. After the Takagi diagonalization of the neutrahno mass matrix at 
M[ = we find two almost massless eigenstates (dominated by the MSSM bino and 
wino components), a pair of nearly degenerate, maximally mixed MSSM higgsinos at 
~ 600 GeV (the third and fourth states) and a pair of nearly degenerate, maximally 
mixed singlino/bino' neutralinos at 949 GeV (the fifth and sixth). The LEP limit on 
the hghtest chargino mass therefore enforces M[ > 55 GeV. 

For understanding the neutralino mixing pattern as a function of M[ it is important 
to remember that the lighter of the two singlino/bino' and the lighter of the two higgsino 
states belong to negative eigenvalues, while the other states to positive eigenvalues. 
When the M[ parameter is switched on, the mixing pattern is more rich since not only 
the singhno/bino', but also the bino and wino states vary considerably, see Fig. 10. As 
a result there are more cross-over zones where mixing is important. In the cross-over 



28 




Figure 10: The neutralino mass spectrum as a function of M[ in the unified scenario 
M[ = Ml. The effective parameter is set to 600 GeV, other parameters as in the 
previous subsection. Right panel is a magnified part of the left one. 

zone around M[ ~ 270 GeV the wino mixes with the heavier higgsino, around 500 
GeV the bino mixes with the heavier higgsino, around 550 GeV the wino mixes with 
the heavier singhno/bino' and in the last zone around 900 GeV the hghter higgsino 
mixes with the hghter singhno/bino' state. This is illustrated in Fig. 10. As the lines 
develop from M[ = GeV, the dominant component of the corresponding state changes 
its nature. For example, along the green line the state starts at M[ = GeV as a 
heavier higgsino, then gradually becomes a wino-dominated (for M( ~ 400 — 600 GeV) 
and finally (for M[ > 600 GeV) a bino'-dominated neutralino. The LSP mass, as we 
increase M[, first increases, then levels off at Mlsp ~ 600 GeV and then decreases 
along with M[. Its nature also changes. It starts as a bino, at M[ ~ 600 GeV gradually 
changes to a higgsino-dominated state and at M[ ~ 800 GeV discontinuously jumps to a 
singlino/bino'-dominated state. For higher values of M[ the lightest neutralino becomes 
mostly singlino. 

5.3.2 Relic density 

In Fig (11) we show the relic density calculation for coupled gaugino masses and /j, — 
600 GeV. In this case the relic density phenomenology is significantly more complex 
than previously. First of all, note that below M[ =0.75 TeV the LSP is predominantly 
bino, with non-zero admixtures from all other states. Above M[ = 0.75 TeV the LSP 
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Figure 11: The relic density across varying M[ with M[ = Mi = M2/2, for fi = 600 GeV, 
rriA = 500 GeV and tan/5 = 5. 

is predominantly singlino with substantial admixtures of bino' and higgsino. Around 
M[ = 0.75 TeV the LSP is predominantly higgsino with a large admixture of both 
singlino and bino. 

If we initially ignore the resonances we can see a general trend in the relic density from 
a large value at low M[, down to a lower value at around M[ = 0.75 TeV and then back 
to larger values at high M[. This is to be expected as this follows the evolution of the 
LSP from bino (that generally gives flcDuh'^ ^ ^wmap) through higgsino (generally 
VLcDAih'^ < ^WMAp) to singlino {QcDMh"^ > ^wmap)- 

Beyond this general structure there are a number of interesting features. Note that 
as M[ increases the LSP mass first increases reaching a maximum of ~ 560 GeV at 
M[ ~ 800 GeV and then falls down crossing all possible s-channel resonances twice. 
Starting from M[ = we first arrive at a little dip in the relic density around M[ = 250 
GeV which is due to the s-channel H2/A resonance. The next resonance due to Z2/H^ 
around M[ = 500 GeV produces only a little wiggle since the LSP has not yet developed 
an appreciable singlino component. The first appreciable dip in the relic density occurs 
around M[ = 0.8 TeV where VLcoNih^ drops to ~ 0.02. Here the LSP has a strong 
higgsino component which enhances the annihilation via the s-channel Z2/H^ resonances 
considerably. Increasing M[ further, the LSP mass increases, going off-resonance (hence 
local maximum in the relic density), until it reaches its maximum of ~ 590 GeV at 
M[ ~ 800 GeV. From now on the LSP mass decreases and its nature becomes singlino- 
dominated. Around M[ = 1.5 TeV it once again hits the Z2/H3 resonance. However, 
this time the LSP is predominantly singlino. Although pure singlino neutralinos do not 
couple to the singlet Higgs, so the resonance is subdominant, they couple strongly to 
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the Z' and annihilate very efficiently. As a result, the rehc density drops to ~ 2 x lO"^. 

The next feature of interest is the kink at M{ = 2.5 TeV. This is where the LSP mass 
drops below threshold for production of HxA in the final state. This backs up our expec- 
tation that annihilation to heavier Higgs states significantly increases the annihilation 
rate of a singlino LSP. 

Prom this point on the relic density profile shows the same essential features as in 
Scenario A. We find a pscudoscalar Higgs resonance at Af{ = 3.5 TeV, the top threshold 
at M[ = 5 TeV, the light Higgs threshold at M[ = 11 TcV, the light Higgs resonance at 
M[ = 14 TeV and the Z resonance at M[ = 19 TeV. The one important difference that 
is worth noting is that in this figure the light Higgs resonance does lower the relic density 
to a point where it agrees with the WMAP-5 measurements. This is due to the doubling 
of A between the two cases. This strengthens the couphng of the singlino-higgsino-higgs 
vertex. 

In our study of Scenario B we can clearly see the effects of increasing the size of yU. 
We can have a heavier singlino which can annihilate to a wider range of final states. 
The singhno also has stronger couphngs to the other Higgs and higgsino states, further 
reducing the relic density. However we see once again that we need to tune the mass of 
the singlino through M[ to fit the relic density, either through a precise balance of the 
singlino/higgsino mixture, or through a careful balance of the singlino mass against the 
mass of a boson that mediates annihilation in the s-channel. 

5.3.3 Direct detection 

In Fig. 12 the spin-independent as well as spin-dependent elastic cross section of the 
hghtest neutralino on proton is shown as a function of M[ . We restrict the range of M[ 
to 0-3 TeV, as beyond this range the cross section falls monotonically. 

Referring to Fig. (10), it is easy to understand the M[ behavior of the cross section. 
For small M[ the lightest neutralino (up to M[ ~ 0.3 TeV) is almost a pure MSSM bino 
and its couplings are roughly M(-independent. As M[ approaches 500 GeV, the LSP 
receives an appreciable admixture of both higgsinos. As a result both spin-independent 
and spin-dependent cross sections rise. However, the spin-dependent cross section being 
sensitive to the combination C34 develops a dip around M[ = 800 GcV until the disconti- 
nuity where two lightest states cross. Above 800 GeV the steady increase of the singlino 
component in the LSP makes the behavior of the cross section resemble the one in the 
previous scenario (for M[ > 2.5 TeV). 

6 Summary and conclusions 

In this paper we have provided an up to date and comprehensive analysis of neutralino 
dark matter within the USSM which contains, in addition to the MSSM states, also 
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Figure 12: The elastic spin-independent (left) and spin- dependent (right) LSP-proton 
cross section as a function of M[ in the unified scenario M[ = Mi . 



one additional singlet Higgs plus an extra Z', together with their superpartners the 
singlino and bino'. We have seen that the extra states of the USSM can significantly 
modify the nature and properties of neutralino dark matter relative to that of the 
MSSM and NMSSM. Using the LanHEP package, we have derived all the new Feynman 
rules relevant for the dark matter calculations. We have also provided a complete 
qualitative discussion of the new annihilation channels relevant for the calculation of the 
cold dark matter rehc density for the neutralino LSP in the USSM. We also discussed 
the elastic scattering cross section for the neutralino LSP in the USSM, including both 
spin-independent and spin-dependent parts of the cross sections, relevant for the direct 
dark matter search experiments. 

We then surveyed the parameter space of the USSM, and discussed quantitatively 
how the nature and composition of the neutralino LSP can be significantly altered com- 
pared to that in the MSSM due to the extra singlino and bino' states, for different 
ranges of parameters. We have considered two approaches to the parameter space: (a) 
holding the MSSM higgsino and gaugino mass parameters fixed, while the mass of the 
extra U(l) gaugino taken free (to complement the collider phenomenology discussed in 
Ref. [18]); (b) the scenario of unified gaugino masses. The Feynman rules were then im- 
plemented into the micrOMEGAs package in order to calculate the relic density for the 
corresponding regions of parameter space. This provides a full calculation of the annihi- 
lation channels including co-annihilation and careful treatment of resonances as well as 
accurately calculating the relic density for an arbitrary admixture of states. In this way 
we extended the analysis of USSM dark matter annihilation beyond the specific cases 
previously studied in the literature. We also performed an equally general calculation 
of the direct detection cross-sections for USSM dark matter for elastic neutralino-nuclei 
scattering. 
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The results show that there are many cases where successful relic abundances may 
be achieved, and in novel ways compared to the MSSM or NMSSM (see for example the 
low mass region in Fig. 11 for M[ < 5 TeV). In general our results also show that the 
inclusion of the bino' state, as well as the lack of a cubic interaction term S^', results 
in a significant change in the dark matter phenomenology of the USSM as compared 
to that of MSSM or NMSSM. Also the neutralino mass spectrum in the USSM may be 
very different from that of the NMSSM as the singlino mass is determined indirectly 
by a mini-see-saw mechanism involving the bino' soft mass parameter M[ rather than 
through a diagonal mass term arising from the cubic S^. The lack of a cubic interaction 
term also restricts the annihilation modes of the singlino, making it dominantly reliant 
on annihilations involving non-singlet Higgs bosons and higgsinos. As the USSM has 
a different Higgs spectrum to the NMSSM, notably in the pseudoscalar Higgs sector, 
the Higgs dominated annihilation channels of the USSM singlino are significantly mod- 
ified with respect to the NMSSM singlino. As Higgs exchange diagrams dominate the 
direct detection phenomenology, the difference in the Higgs spectrum and the singlino 
interactions results in significant differences in the direct detection predictions as well. 

In conclusion, the USSM, despite its modest additional particle content compared 
to the MSSM or NMSSM, leads to a surprisingly rich and interesting dark matter phe- 
nomenology which distinguishes it from these models. The other states which are nec- 
essary in order to make the model anomaly free, and which we have neglected here, 
can only add to the richness of the resulting phenomenology, but the qualitatively new 
features that we have found in the USSM will remain in any more complete model. 
Nevertheless it would be interesting to study the effect of the additional states present, 
for example, in the EgSSM in a future study. 
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A Higgs boson masses 

In general the neutral CP-even Higgs fields h, H, S mix. The mass matrix takes the 
form (see the first paper in [27]) 

Mf^ Mf^ Mfg 

MLn = I Ml Ml Ml I (A.l) 



Ml Ml Ml 



where 



Ml = ^^;2sin2 2/3+^^^^;2cos^2/3 + ^;V(giCOS^/3 + g2sin2/3)2 + An, 

Ml = Ml={^-3-^^y'^r.Ap+^lv\Q,-Q,)x 

X (Qi cos^ /3 + Q2 sin^ /3) sin 2/3 + A12 , 
Ml = ^(^^5+f^^-^V'«in'2/3+£(Q2-Qi)Vsin2 2/3 + A^22, 



sin2/3 " \^ 4 2 y ^4 

^13 = Ml = - sin 2(5 + A^t; + g'^ {Qi cos^ /3 + Q2 sin^ /3) Qst; + A13 , 
v2 

^la = M|2 = --^t;cos2/3+^(Q2-gi)g5^^;5sin2/3 + A23, 

= ^^^v' sm2p + g['Qlvl + Aess (A.2) 
where the one loop-corrections Aeij are expressed as 

Ae22 = A22 + Aa-A^ (A.3) 

As 

Ae33 — A33 

Vs 

Ac 

Aea = Aa-A^- — + A3 

Vs 

in terms of Aij, As, Ap, Aa and A3 given explicitly in Ref. [33] (note that the expression 

m| mi 

for X in this paper should read K — F — ^ log — hr^)- 



B Feynman rules 

All the Feynman rules presented here are given in terms of the interaction of mass 
eigenstates. As a result the Feynman rules reference many matrices that rotate from 
the interaction eigenstate to the mass eigenstate basis. We briefly summarise them here 
for ease of reference: 
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Figure B.l: The ^fx?Xn vertex given in Eq. (B.l) 

• a, - Z,Z' mixing matrix that transforms from the Z' eigenstates to the Z\ 2 
mass eigenstates, as defined in Eq. (10) 

• - Higgs mixing matrix from the mass eigenstate basis to the the interaction 
eigenstate basis, defined in Eq. (18). 

• Nij - neutrahno mixing matrix, defined in Eq. (26) 

• t/jj, Vij - standard chargino mixing matrices as in the MSSM [4]. 

• C/J - squark or slepton mixing matrix, defined in Eq. (31). 

Feynman rule for the vertex shown in Fig. B.l: 



A COS (7^y 
Z COS uw 



(B.l) 



Feynman rules for the Xk'Xd^^ vertex shown in Fig. B.2 

ig,r KkPL + C^Pr) (B.2) 
tg,r {C^^Pl + ClCPn) (B.3) 
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(a) 




W-'" (b) 




Figure B.2: The XfeX?^^^ S^^®" (^) (^-2) and (b) Eq. (B.3) 




Figure B.3: The x°Xn-^fc vertex given in Eq. (B.6) 



where 



^1 



V2 

1 

3^ 



Feynman rule for the x'lXn-^k vertex shown in Fig. B.6: 
where 



Rni — 



nl 



T, 



nl 



+(Z ^ n) 

^(iV„2 - tan^w'iV„i)iV;4 - g[Q2Nn(iNi^ + ^^nsiV^s 
+(Z ^ n) 

-g'lQsNn&Ni^ + :^^n3A^/4 + ^ n) 



(B.4) 
(B.5) 

(B.6) 



(B.7) 

(B.8) 
(B.9) 



36 



Figure B.4: The XiXn^ vertex given in Eq. (B.IO) 



(a) 




-H- 



(b) 




- -H+ 



Figure B.5: The X^X?//^ vertex given in (a) Eq. (B.14) and (b) Eq. (B.15) 



Feynman rule for the x?Xn^ vertex shown in Fig. B.IO: 

[(i?;; sin P + S'l* cos P) cos + r;* sin 0] Pl 
- [{R'^i sin /3 + S'^i cos (3) cos + T^, sin 0]P« 



where 



nl 



nl 



+ (/ ^ n) 

A 



Feynman rules for x^ X?-^^ shown in Fig. B.5: 

37 



(B.IO) 



(B.ll) 

(B.12) 
(B.13) 

(B.14) 



Figure B.6: The ^^^X? vertex given in (a) Eq. (B.18) and (b) Eq. (B.19) 



where 



R 



,"L 



V* 



Wi + ^W2 + A^ntanM 



Hk = ^2COS/3 

+g[V2 COS (3Q2N:^V:^ + A sin/5iV* 
i?;'/ = ^2sin/5 Ni:,Uui-^{Ni2-Nat^nt 

+g[\f2sml3QiNiQUk2 + Xcos (3Ni^Uk2 
Feynman rules for the qqs)(^ vertex shown in Fig. B.6 

I [GI^Pl + GlfPn] 
For up- type quarks: 



G:t = -^/2 



92 i\N:^ + i tan OwN^^ ) + g'^QqN*^ 



(B.15) 



U- - 

V2Mw sin P 



G7 = ^ 



g2- tan OwNn - g[QuNi6 



jj2s _ 92mu, 

V2Mw sin p 



(B.16) 

(B.17) 

(B.18) 
(B.19) 



N* TT"^^ 



(B.20) 



For down- type quarks: 
Gt? = -V2 



92 QiV* - ^tan^i^TV*^ - g'^Q^N^, 



* V2Mwcos(3 * 



tan ewNa + g[QaNie 



Uf - ^^1. 

* V2Mw cos P * 



(B.21) 
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